Oscillations of the large-scale circulation in turbulent Rayleigh-Benard convection: 
the off-center mode and its relationship with the torsional mode 
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We report an experimental study of the large-scale circulation (LSC) in a turbulent Rayleigh- 
Benard convection cell with aspect ratio unity. The temperature-extremum-extraction (TEE) 
method for obtaining the dynamic information of the LSC is presented. With this method, the 
azimuthal angular positions of the hot ascending and cold descending flows along the sidewall are 
identified from the measured instantaneous azimuthal temperature profile. The motion of the LSC 
is then decomposed into two different modes based on these two angles: the azimuthal mode and 
the translational or off-center mode that is perpendicular to the vertical circulation plane of the 
LSC. Comparing to the previous sinusoidal-fitting (SF) method, it is found that both the TEE and 
the SF methods give the same information about the azimuthal motion of the LSC, but the TEE 
method in addition can provide information about the off-center motion of the LSC. The off-center 
motion is found to oscillate time-periodically around the cell's central vertical axis with an ampli- 
tude being nearly independent of the turbulent intensity. It is further found that the azimuthal 
angular positions of the hot ascending flow near the bottom plate and the cold descending flow near 
the top plate oscillate periodically out of phase by 7r, leading to the torsional mode of the LSC. 
These oscillations are then propagated vertically along the sidewall by the hot ascending and cold 
descending fluids. When they reach the mid-height plane, the azimuthal positions of the hottest 
and coldest fluids again oscillate out of phase by ir. It is this out-of-phase horizontal positional 
oscillation of the hottest and coldest fluids at the same horizontal plane that produces the off-center 
oscillation of the LSC. A direct velocity measurement further confirms the existence of the bulk 
off-center mode of the flow field near cell center. 

The paper has been submitted to J. Fluid Mech. 



I. INTRODUCTION 

The convection phenomenon is ubiquitous in nature 
and in many engineering applications. A simple but 
paradigmatic model that has been widely used to study 
the convection phenomenon for more than a century is 
the turbulent Rayleigh-Benard (RB) convection, which 
is a fluid layer heated from below and cooled on the top 
(for recent reviews, see @, Hl|). The dynamics of the 
RB system is determined by its geometry and two di- 
mensionless control parameters: The Rayleigh number 
Ra = ftgH 3 AT ' jvn and the Prandtl number Pr = v j 'k, 
where g is the gravitational acceleration, H the height of 
the cell, AT the temperature difference across the cell, 
and (3, v, and k, respectively, the thermal expansion coef- 
ficient, the kinematic viscosity, and the thermal diffusiv- 
ity of the fluid. The geometry of the system is described 
by its symmetry and the aspect ratio of the convection 
cell F. For a cylindrical cell, r = D/H with D being the 
cell's diameter. 

A prominent feature of turbulent RB system is the 
presence of the large-scale circulation (LSC), which is 
self-organized from thermal plumes that erupt from the 
top and bottom thermal boundary layers [28] . We call it 
"large-scale" because it is a single cellular structure that 
spans the height of the cell, at least in cells with aspect 
ratios close to one. In an axial symmetric configuration, 
the LSC is found to exhibit stochastic azimuthal mean- 
dering and small-probability events of cessations and re- 
versals (see and references therein). In addition, the 



horizontal motions of the LSC near the top and bottom 
plates are found to oscillate periodically out of phase by 
7T. This twisting oscillation of the LSC's vertical circula- 
tion plane is referred to as the torsional mode of the LSC, 
which was first found by Funfschilling & Ahlers 

In addition to the torsional oscillation of the LSC, there 
also exists a well-defined low-frequency oscillation in both 
the temperature and velocity fields, which has been long 
observed in turbulent RB ex per iments using various flu- 
ids (see, e.g., @, 0, [H, H, Ujl HI)- As oscillation is 
a common phenomenon in closed flow systems, under- 
standing the nature and the origin of the low-frequency 
oscillation in the RB system should be of general inter- 
est. To understand this oscillation, Villermaux [13] pro- 
posed the coupled boundary layer model and the two 
main predictions of the model are: (1) Plumes are emit- 
ted periodically and (2) hot and cold plumes are emitted 
alternatively. The two predictions are not independent 
of each other. Physically speaking, thermal plumes can 
emit periodically but not necessarily alternatively. How- 
ever, they cannot emit both alternatively and nonperi- 
odically. This is because the alternative emission of hot 
and cold plumes implies a triggering mechanism of plume 
generation, i.e., hot plumes cause the emission of cold 
ones and vice versa. Therefore, the emission of one type 
of plumes should wait for the arrival of the other type. 
Since there exists a typical timescale for thermal plumes 
to travel from one plate to the other, the alternative 
emission process should quickly lead to a synchronized 
and periodic emission of hot and cold plumes from the 
bottom and top plates respectively. Some experiments 
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based on single-point or two-dimensional (2D) measure- 
ments of the temperature and velocity fields appear to 
support Villermaux's model [1, [H, Ell, Hif- Some 
other works, however, pointed out that periodic plume 
emission is not necessary for the periodicity of the sys- 
tem [13, Ell- Recently, Ahlers et al. Q conjectured that 
the low-frequency oscillation of the system presumably 
is due to the torsional oscillation of the LSC. However, 
how the torsional mode generates the temperature and 
velocity oscillations at the mid-height plane of the cell is 
unclear, since, without the off-center oscillation, a simple 
twisting oscillation near the top and bottom plates would 
cancel out at the mid-height plane due to symmetry. 



In a recent experimental study of the three- 
dimensional (3D) structure of temperature oscillations, 
Xi et al. f32| have presented convincing evidences that 
hot and cold thermal plumes are emitted neither period- 
ically nor alternatively, but continuously and randomly, 
from the top and bottom plates. Thus invalidate the 
two predictions of the Villermux's model. They further 
showed that the low-frequency oscillation at the cell's 
mid-height plane is due to the off-center oscillation of the 
LSC. It should be mentioned that a previous single-point 
velocity measurement has shown that at the cell center 
the strongest velocity oscillation is along the direction 
perpendicular to the LSC plane and that the strength of 
this oscillation decays away from the cell center towards 
the plates PU . More recently, a study of the 3D spatial 
structure of the velocity field has shown that the veloci- 
ties along the axis perpendicular to the LSC plane and at 
the cell's mid-height plane correlate strongly with each 
other and have a common phase across the cell's entire 
diameter (24|. Both these results imply the existence of 
the off-center motion. However, the nature of this mo- 
tion and its relationship with the torsional oscillation of 
the LSC have not been revealed, which are among the 
objectives of this paper. 



The rest of the paper is organized as follows. We de- 
scribe the experimental setup and conditions in Section 
2.1. In Section 2.2, we describe in detail a method for 
extracting the azimuthal angular positions of the hottest 
and coldest fluids along the sidewall and the LSC's cen- 
tral line from the measured instantaneous azimuthal tem- 
perature profile and provide a validation of this method. 
Comparisons with the previous sinusoidal-fitting method 
will also be made. The experiment results are presented 
and analyzed in Section 3, which is divided into three 
parts. In Section 3.1 we present a detailed study of the 
off-center oscillation of the LSC. Section 3.2 discusses 
the relationship between the off-center mode and the tor- 
sional mode of the LSC and Section 3.3 presents results 
from a direct velocity measurement of the bulk off-center 
mode of the flow field. We summarize our findings and 
conclude in Section 4. 
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FIG. 1: (a) Top- view of the convection cell and coordinate sys- 
tem used. Also shown is an actual instantaneous azimuthal 
positions ah and a c of the hot ascending and cold descending 
flows, corresponding to the temperature profile T(a) shown 
in (c). The z-axis (not shown in the figure) is along the cell's 
central vertical axis. The numbers 1, 2, . . . , 8 show the loca- 
tions of the sidewall thermistors, (b) Side-view of the cell with 
measured circulation path of the LSC reproduced from Ref. 
[23I ]. The three dashed lines show the top, middle and bottom 
planes where the thermistors are placed, (c) An example of an 
actual instantaneous temperature profile (circles) measured 
by the 8 thermistors at mid- height plane at Ra — 5.5 x 10 9 . 
The solid line represents a quadratic fit to the highest tem- 
perature reading and the two temperature readings adjacent 
to it. The dashed line represents a quadratic fit to the lowest 
temperature reading and the two temperature readings adja- 
cent to it. The two vertical arrows indicate the fitted peak 
positions of aj, and a c , respectively. The dotted line shows a 
sinusoidal fit to the profile. It is seen from the figure that, in 
this case, the sinusoidal fit to the profile shifts aj, to the right 
(anticlockwise in (a)) and shifts a c to the left (clockwise in 
(a)). 



II. EXPERIMENTAL SETUP, METHODS AND 
PARAMETERS 

A. The convection cell and experiment conditions 

The experiment was carried out in a cylindrical cell 
with its top and bottom plates made of 1.5-cm-thick cop- 
per, the sidewall of a 5-mm-thick Plexiglas tube and wa- 
ter as convection fluid (30j . The inner diameter D and 
the height H of the cell are both 19.0 cm and hence its as- 
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pect ratio V is unity. Twenty-four thermistors (Omega, 
44031) with a diameter of 2.5 mm and an accuracy of 
0.01° were placed in blind holes drilled horizontally from 
the outside into the sidewall with a distance of 0.7 mm 
from the fluid-contact surface. These thermistors are dis- 
tributed in three horizontal rows at distances H/A, H/2, 
and 3H/4 from the bottom plate, which are denoted as 
the bottom, middle and top planes (the three dashed lines 
in figure QJb)), respectively, and in eight vertical columns 
equally spaced azimuthally around the cylinder (figure 
[lja)). A multichannel multimeter was used to record 
the resistances of the 24 thermistors at a sampling rate 
of 0.29 Hz, which are converted into temperatures using 
calibration curves. Then the azimuthal temperature pro- 
file T(a) at the three heights can be obtained, where a 
is the azimuthal angle referenced to location 1. During 
the measurements, the mean temperature To of the bulk 
fluid was kept at 31°C and hence Pr = 5.3. The measure- 
ments covered six values of Ra, ranging from 9.0 x 10 8 
to 6.0 x 10 9 , and lasted 70 to 750 hours. 

In Ref. [32| the cell was tilted by 2° to study the ori- 
gin of the temperature oscillations. This is because in 
addition to the twisting and off-center sloshing motion, 
the orientation of the LSC also meanders randomly. By 
tilting the cell and thus locking the LSC orientation, one 
can remove the stochastic meandering from the signal 
and separate the complicated phenomena produced by 
the different types of motions. This enabled us to study 
the phase relationships between temperature oscillations 
at various locations. In this work, our focus is on the 
relationship between the twisting motion near the top 
and bottom plates and the off-center motion at the mid- 
height plane. As the LSC meanders azimuthally as a 
whole across the height of the cell [23] j the phase rela- 
tionship between these two types of motions can be stud- 
ied even with the azimuthal meandering present. There- 
fore, unless stated otherwise, all measurements in the 
present work were made with the cell levelled (to within 
0±0.06°). 

As stated above, the temperature measurements in the 
present case are made with thermistors embedded in the 
sidewall, whereas in Ref. [Hj] the off-center motion was 
obtained from thermistors placed in fluid. It is known 
that the sidewall acts as a low-pass filter, thus the ther- 
mistors embedded in the sidewall are not sensitive to the 
high-frequency signals. Spatially, they actually sense the 
integrated signal over a finite area of the sidewall, which 
leads to a lowered strength of the off-center oscillation 
as compared from that using the in- fluid probes. Never- 
theless, as we will see below, the basic pictures obtained 
from the two cases are essentially the same. As one of the 
objectives of the present work is to understand the rela- 
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tionship between the off-center oscillation and the twist- 
ing oscillations near the plates, we use data measured 
simultaneously by the 24 in-wall probes from the three 
heights (there are only 8 in-fluid probes and that mea- 
surements can only be made at one height at a time). 
Thus, unless stated otherwise, all results presented in 
this paper were obtained from the in-wall probes. 



B. The temperature-extremum-extraction method 

An often-used method based on the multi-probe tech- 
nique for extracting the dynamic information about the 
LSC motion is to fit the temperature azimuthal profile 
using a sinusoidal function, i.e. Tk = T a + A'cos(kir/4 — 
4>'),k = 0,1,. ..,7, where T a is the azimuthal average 
of the 8 temperature readings, A' is a measure of the 
strength of the LSC and d/ is the LSC's orientation 
(see, e.g. d H i, [ill, US S3])- This sinusoidal-fitting 
(SF) method has been very successful in the study of 
the azimuthal motion of the LSC, including rotations, 
cessations and reversals (see [3] and references therein). 
However, the SF method requires the separation between 
the hottest and coldest azimuthal positions to be it, i.e., 
the obtained LSC's central line is forced to always pass 
through the cell's central vertical axis, which, as we shall 
show below, is not always the case. 

Here, we introduce the temperature-extremum- 
extraction (TEE) method that determines first the 
hottest and coldest azimuthal positions of the bulk fluid 
and then the central line of the LSC band. The TEE 
method has been described very briefly in Ref. [32] ]. 
In this section, we give its detailed description, valida- 
tion and comparison with the SF method. To illustrate 
this mehtod, figure UX C ) shows an example of the actual 
instantaneous temperature readings (circles) from the 8 
thermistors at mid-height plane. The hottest (coldest) 
azimuthal position of the bulk fluid along the sidewall 
was determined by making a quadratic fit around the 
highest (lowest) temperature reading. In practice, this 
requires only the local maximum (minimum) tempera- 
ture and the two temperatures directly adjacent to it. 
This is because a quadratic function has 3 degrees of 
freedom and thus can be uniquely determined by only 3 
data points. We label these azimuthal positions ai, a 2 
and c*3 with a 2 being the position of the local maximum 
(minimum) temperature. The position a.h (a c ) at any 
plane can then be found by solving analytically the 3 
quadratic equations between T(afe) and a^ (k — 1, 2, 3), 
i.e., 
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FIG. 2: The averaged sidewall-temperature-profile [T(a) — 
To] normalized by the temperature difference AT around the 
fitted (a) hottest and (b) coldest azimuthal positions obtained 
at Ra — 5.5 x 10 9 . Solid lines are quadratic functions. 



In this paper, the subscripts 'h' and 'c' are used to denote 
those quantities associated with the hot ascending and 
cold descending flows, respectively, and the subscripts 'i', 
'm' and '6' for the top, middle and bottom planes, respec- 
tively. To test the validity of the quadratic fit, figure [5] 
shows the normalized sidcwall-tcmperaturc-profile. Each 
point is an average of [T(ctk) — To] /AT (k — a±, a2, a^) 
in a bin with a small range around a — ah. m (figure 
[DJa)) or a — a c . m (figure [^b)). The data around the 
hottest (figure HJa)) and coldest (figure HJb)) positions 
are both in good agreement with the quadratic functions 
(solid lines), indicating that the quadratic function is in- 
deed a good representation for the temperature distribu- 
tion around the hot ascending and cold descending flows. 
With this method, the three rows of thermistors can thus 
provide simultaneously the azimuthal positions of the hot 
ascending and cold descending flows at the three heights, 
which are denoted as au.ti cth,m, and oth,b, and a c .t, a c . m , 
and a Ci b. With the obtained hottest and coldest posi- 
tions, the line connecting the two positions is the central 
line of the LSC band. The orientation of this line is the 
orientation of the LSC, the distance between this line and 
the cell's central vertical axis is defined as the off-center 
distance (see figure QJa)) and A — (T(ah) — T(a c ))/2 is 
used to characterize the strength of the LSC, where the 
factor 1/2 is used in order to compare the obtained A 
with that obtained from the SF method. We denote <pt, 
4> m , and 4>b, A t , A m , and Ab, and d t , d m , and db as the 
orientations, the strengths and the off-center distances of 
the LSC at the top, middle, and bottom planes, respec- 
tively. Therefore, based on the TEE method, the motion 
of the LSC can be decomposed into two different modes: 
The azimuthal mode and the off-center mode. 

We now compare the orientations and the strengths of 
the LSC obtained using the TEE method and those ob- 
tained using the SF method. Figure O^a) shows the time 
traces of 4>m and 4>' m , which are obtained using the same 
data set measured by the embedded thermistors. One 
sees that the time traces are very similar to each other. 
This similarity can be characterized by the coefficient of 
the cross-correlation between the two quantities. The 
cross-correlation function between two variables a and b 
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FIG. 3: (a) Time traces of the orientation 4> m of the LSC cen- 
tral line obtained using the TEE method and the orientation 
4>'m obtained from the SF method, (b) The cross-correlation 
functions between the orientations obtained using the two dif- 
ferent methods for the top (solid line), middle (dashed line) 
and bottom (dotted line) planes, (c) PDF of the two orienta- 
tions 4>m and 4>' m . All data were measured at Ra = 5.5 x 10 9 . 

is defined as 

C a , b (r) = {(a(t + t) - {a)){b{t) - {b))) / a rms b rms , (2) 

where a rms and b rms are standard deviations of a and 
&, respectively, and (■ • ■ ) represents the temporal aver- 
age. When a = b, C a , ai denoted as C Q , is the auto- 
correlation function of the variable a. Figure [3]Jb) shows 
the cross-correlation functions between <p and (/)' at the 
three heights. It shows that all these three functions have 
a positive and strong correlation with the coefficient be- 
ing all larger than 0.85 at time lag r = 0. (The subpeaks 
°f C^t.^i ano - ^<Pb,tl>' b correspond to the twisting oscilla- 
tions of the LSC near the top and bottom plates.) It is 
further found that this strong positive correlation exists 
for all values of Ra investigated. Figure [3£c) shows the 
probability density functions (PDF) of 4> m and <j)' m . One 
sees that the two PDFs essentially collapse on top of each 
other except that the distribution of 4>' m seems to be a 
little smoother than that of (f> m . 

Figure Ufa) shows the time traces of A m and A' m . It 
is seen that the two time traces are again very similar to 
each other except that A m is on average ~ 5% larger than 
A' rn . Figure BJb) shows the cross-correlation functions 
between A and A' at the three heights. It similarly shows 
that all three functions have a strong positive correlation 
with their peaks all located at time lag r = and the 
cross-correlation coefficient C(r = 0) all being nearly 0.9. 

The above results imply that, as far as the orienta- 
tion and the strength of the LSC are concerned, both 
the SF and TEE methods give the same results, but the 
TEE method in addition determines the LSC's off-center 
motion, which is missed by the SF method. This is il- 
lustrated in figures [5]Ja) and (b). It is also clear that the 
LSC's orientation is determined mainly by the hottest 
and coldest azimuthal positions of the bulk fluid along 
the sidewalk This is not surprising, since one can see 
from figure [TJc) that the sinusoidal fit to the temper- 
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FIG. 4: (a) Time traces of the normalized strength A m /AT 
of the LSC obtained using the TEE method and the normal- 
ized strength A' m / AT obtained from the SF method, (b) The 
cross-correlation functions between the strengths obtained us- 
ing the two different methods for the top (solid line), middle 
(dashed line) and bottom (dotted line) planes. The data were 
measured at Ra = 5.5 x 10 . 
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FIG. 5: (a) Schematic diagram showing the definitions of the 
orientation (j> and the off-center distance d of the LSC obtained 
using the TEE method, (b) Schematic diagram showing the 
definition of the orientation <j> of the LSC obtained using the 
SF method. In this method, the off-center distance is forced 
to be always zero. 
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FIG. 6: (a, b) Time traces of the normalized off-center dis- 
tance d m /D from thermistors (a) embedded into the sidewall 
and (b) placed in fluid, (c) PDF of d m /D from the in- wall 
(solid circles) and in-fluid (open circles) thermistors. The 
data of the in-fluid thermistors are taken at Ra — 5.0 x 10 9 
with a sampling rate of 0.76 Hz and that of the in-wall probes 
at Ra = 5.5 x 10 9 with a sampling rate of 0.29 Hz. 



ature profile shifts the hottest azimuthal position oth,m 
right or anticlockwise and shifts the coldest azimuthal 
position a c . m left or clockwise and the changes of the 
obtained orientation of the LSC due to the two shifts 
roughly cancel each other. 



III. RESULTS AND DISCUSSION 

A. The off-center oscillation of the LSC at 
mid-height plane 

The off-center distance d m of the LSC's central line 
at the mid-height plane is used to study the off-center 
motion of the LSC at that plane. Figure G^a) shows a 



typical time trace of the measured d m normalized by the 
cell's diameter D. It is seen from the figure that d m /D 
fluctuates around 0. We further found that the temporal 
averages of d m /D for all six values of Ra are nearly 0. In 
Ref. 32[ a similar trace measured with the in-fluid ther- 
mistors in a tilted cell has been shown. For comparison, 
we plot in figure [^b) a time trace of d m /D measured 
under the same conditions as in Ref. [32j but with the 
cell leveled (0 ± 0.06°). It shows that d m /D measured 
from the in-fluid probes has a stronger strength and its 
oscillation looks more periodic when compared with that 
obtained from the in-wall thermistors. Nevertheless, the 
physical pictures revealed by the two time traces are the 
same, i.e., the LSC's central line oscillates horizontally 
around the cell's central vertical axis. Figure [^c) shows 
the PDF of d m /D obtained from the in- wall thermistors 
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This result further indicates that the previously-observed 
temperature oscillations near the mid-height of the cell 
indeed originates from the off-center oscillation of the 
LSC. Taken together, the physical picture behind figure 
[6]is that the LSC's central line oscillates time-periodically 
around the cell's central vertical axis with an amplitude 
nearly independent of the turbulent intensity. 
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FIG. 7: (a) The i?a-dependence of the standard deviation of 
dm/D from the in- wall (solid circles) and in- fluid (open cir- 
cles) thermistors. Two dashed lines show the mean value of 
d m / D averaged over all Ra investigated, (b) Normalized oscil- 
lation frequency H 2 fo/K. obtained from the frequency power 
spectra of d m as a function of Ra. The solid line is the power- 
law fit. 



(solid circles). The PDF is found to be flatter than a 
Gaussian function, i.e. the flatness is 2.3. The PDFs ob- 
tained at other Ra share the same features as that shown 
in hgurelHfc), except for the two lowest Ra which is prob- 
ably due to the limited temperature contrast along the 
sidewall. For comparison, we also plot the PDF of d m /D 
obtained from the in-fluid thermistors (open circles), ft 
is seen that the PDF from the in-fluid probes is even flat- 
ter (flatness 1.9) than that from the in-wall thermistors, 
reflecting a larger range of oscillation. 

As d m /D oscillates around its mean value of 0, its 
standard deviation (d m ) rms /D can be used as a mea- 
sure of the amplitude of this oscillation. Figure [Ha) 
shows the i?a-dependence of (d m ) rms / D obtained from 
the in-wall (solid circles) and in-fluid (open circles) ther- 
mistors. One sees that (d m ) rms /D depends weakly on Ra 
for the present range of Ra with (d m ) rms / D from the in- 
fluid probes being ~ 60% larger than that obtained from 
the in-wall thermistors. In Ref. [32| it has been shown 
that the off-center or sloshing motion of the LSC exhibits 
a well-defined time-periodic oscillation. The prominent 
peak near /o in the frequency power spectra of d m (see 
the second curve from bottom in figureJSJa)) corresponds 
to this periodic oscillation. In figure [Tfb) we plot the 
i?a-dependence of the normalized frequency correspond- 
ing to this prominent peak obtained from the spectra of 
d m . The solid line in the figure represents the power-law 
fit to the data: H 2 f /n = 0.12Ra OA9±Q 02 , with the scal- 
ing exponent being in good agreement with the previous 
temperature measurements in water [El. ITU. 17 ,_ 23|. mer- 
cury [25[ and low-temperature helium gas [J lj 1J, [2(| ■ 



B. The relationship between the off-center mode 
and the torsional mode of the LSC 

Figure |5[a) shows, from top to bottom, the frequency 
power spectra of 4>t, 4> m , 4>b, d tl d m , and db- For the 
orientations of the LSC obtained at the three heights, the 
prominent peak near /o representing periodic oscillation 
can be seen clearly for (j) t and 4>b, but is absent for the 
mid-height plane (4> m \- These are consistent with those 
observed by Ref. IllJ and correspond to the twisting 
mode of the LSC [T(|. We further found that although 
4> m will also exhibit oscillation when the cell is tilted by 
several degrees, the strength of this oscillation obtained 
at the mid-height plane is much weaker than those of <j>t 
and 4>b, i-e., the peak height of the spectra of </> m at /o is 
much smaller than those of <j) t and 4>b- However, for the 
off-center distance, the situation is opposite. In figure 
[5Ja) one sees that the strength of the oscillation peak for 
d m at the mid-height plane is much stronger than those 
of dt and db ■ 

To understand this difference and to find out the rela- 
tionship between the off-center mode at the mid-height 
plane and the torsional mode of the LSC near the top and 
bottom plates, we examine the frequency power spectra 
of the azimuthal positions of the hot ascending and cold 
descending flows at the three heights, as both the ori- 
entation (f> and the off-center distance d of the LSC are 
obtained from these positions. Figure &b) shows that 
the spectra of a c ,t and ah,b exhibit an oscillation peak 
near Jo- Whereas the spectra of ah,t and a c ^ both give 
very faint oscillations, which is due to a tilted ellipse- 
like circulation path of the LSC when viewed from the 
side (Qiu & Tong 2001b; Sun & Xia 2005; Sun et al. 
2005b) . The thermistors embedded into the sidewall can 
not feel accurately the hot ascending flow of the LSC at 
the top-plane and the cold descending flow of the LSC at 
the bottom-plane, since they are far away from the side- 
wall at the respective heights (figure [Ub)). Therefore, at 
the top and bottom planes, the in-wall probes may not 
be able to measure the orientation and the off-center dis- 
tance of the LSC as accurately as those at the mid-height 
plane. This suggests that the relatively weak oscillation 
peaks of dt and db shown in figure EJa) are mainly due 
to the oscillations of a Cj t and cth,b- The surprising thing 
shown in figure [Sfb) is that the spectra of the azimuthal 
positions ah, m and a C;m of the hot ascending and cold 
descending flows at the mid-height plane both exhibit a 
well-defined time-periodic oscillation with the oscillation 
peak located at /o and with the same strength as those 
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FIG. 8: From top to 
bottom: Power spectra 

Of (a) (j>t, <j)m, 4>b, dt, 

d m , and d b , and (b) a h ,t, 
Oth.m, Oth,b, Q!c,t, a c ,m, 
and a c> t- All data have 
been shifted vertically for 
clarity. 
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FIG. 9: Time traces of 

(a) ah,b and a h , m and of 

(b) a c ,t and a c , m - The 
cross-correlation func- 
tions between (c) Oih,m 
and oth.b and between 
(d) a c , m and a c ,t- All 
data were obtained at 
Ra = 5.5 x 10 9 . 



of a.h,b and ct c ,t- This implies that the azimuthal po- 
sitions of the hot ascending and cold descending fluids 
both oscillate periodically along the sidewall irrespective 
of whether they are near the plates or at the mid-height 
of the cell. 

Figure [H£a) shows typical time traces of Uh.b and 
ah,m measured at the same time. It is seen that the 
two time traces are rather similar to each other except 
for a several-second delay between them. The cross- 
correlation function C ah mi a h b shown in figure[5]Jc) quan- 
tifies this similarity. It shows that «h,m and othfi have a 
strong positive correlation with the main peak located at 
t ~ 6 ± 1 s. This positive time-delay indicates that the 
azimuthal motion of ah, m lags that of which is easy 
to understand since the hot ascending flow rises up from 
bottom. Similar situation can be seen for the relation be- 
tween a C; t and a C;m . Figure EJb) shows the time traces 
of a c j and a c ,m, which again are quite similar to each 
other. Figure [U^d) plots C a a t , and it is found that 
otc.m, correlates strongly with a Ci t with a positive time 
delay r ~ 5 ± 1 s indicating that the azimuthal motion of 
a c ,m lags that of a c .t- Here we note that the time delay 



of the main peak of C Qc miQ . c t is essentially the same as 
that of C ah m ,a h b within experimental uncertainty, this is 
due to the fact that the distance between the middle and 
the top planes is the same as the distance between the 
middle and bottom ones. Note also that the subpeaks of 
Ca h m .a h b and C ac mjQc t correspond to the periodic oscil- 
lations shown in figure [5Jb). These results further imply 
that the hot ascending and cold descending flows go up 
and fall down coherently, thus propagate their azimuthal 
positional oscillations along the sidewall. 

The torsional mode of the LSC can be revealed by 
the horizontal motions of the hot ascending and cold de- 
scending flows near the top and bottom plates. Here, we 
define 

Uij = datij/dt, (3) 

as the azimuthal angular velocity of a^j, where i = h 
or c and j — t, to, or b. Figure [TUf a) shows the cross- 
correlation function C Ue tMh b between cj C) j and u>h,b (blue 
line). It is seen that C UJc tMh b oscillates and w c t anticor- 
relates with u>h,b with a strong negative peak located at 
t = 0. The negative peak at r = indicates that a c ,t ro- 
tates clockwise when aih,b rotates anticlockwise and vice 
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© : hot ascending flow (§) : cold descending flow 

FIG. 10: (a) The cross-correlation functions C UettUhb (blue 
line) and C^ am ^ h m from the in- wall (red line) and in-fluid 
thermistors (dark green line). The data from in- wall ther- 
mistors are taken at Ra — 5.5 x 10 9 and that from in-fluid 
probes at Ra = 5.0 x 10 9 . (b) Schematic diagram of the az- 
imuthal positions of the hot ascending and cold descending 
flows when the off-center oscillation is at an extremum. The 
curved arrows show the directions of the angular velocities 
and the horizontal arrow shows the moving direction of the 
LSC's central line, (c) Same as (b) but half an oscillation 
period later. 



versa, implying the twisting motion of the LSC. Recall 
that the azimuthal motions of o.h, m and a C;m follow co- 
herently those of cth.b and a Cj t, respectively, with the 
same time delay. The azimuthal angular velocities of 
ah^m and a Cjm would hence be expected to exhibit the 
same behavior as those of oth,b and a c ,t, and this is indeed 
observed from figure [Tol a), which shows that C Ue m ^ h m 
(red line) is the same as C u)c ttUh b (blue), i.e. periodic os- 
cillations and a strong negative peak at r = 0. Here, we 
also plot C UJc m >u>h m obtained from the in-fluid thermis- 
tors (dark green line). It is seen that C Wc Uh m obtained 
from the in-fluid thermistors exhibits the oscillation with 
much stronger periodicity than that obtained from the 
in- wall thermistors (red line). However, the basic picture 



is the same for both cases. These results indicate that the 
azimuthal positional oscillations of the hot ascending and 
cold descending flows at the mid-height plane are out of 
phase with each other by 7r, leading to the off-center mode 
of the LSC at the mid-height plane. The mechanism of 
how this off-center mode is generated is illustrated in fig- 
ures [TO^b) and (c). Figure [TOT b) shows that a c . m rotates 
clockwise when a.h,m rotates anticlockwise, and the line 
connecting them thus moves to the right. Half an os- 
cillation period later (figure flWc)). the motion of a c irl 
changes to anticlockwise while that of oth^m. to clockwise. 
This makes the LSC central line move to the left. The 
periodic occurrence of this process thus generates the off- 
center oscillation of the LSC and produces the prominent 
peak near fo on the spectra of d m . On the other hand, as 
the orientation <j) m equals to (a.h, m + &c.m)/2 plus a con- 
stant, the anticorrelation between LOh, m and uj c . m would 
cancel out the periodic oscillations, and hence it is not 
surprising to see in figure [5]Ja) that an oscillation peak is 
absent for the spectra of <p m . 



C. Direct velocity measurement of the off-center 
motion 

We present in this subsection direct evidence of the 
off-center oscillation of the LSC from particle image ve- 
locimetry (PIV) measurement of the horizontal velocity 
field at the mid-height plane in a sapphire cell. Both the 
sapphire cell and the horizontal velocity measurement 
using the PIV technique have been described in detail 
previously [3ll . l34j and hence we outline only their main 
features here. The sapphire cell consists of two sapphire 
discs with thickness 5 mm as the top and bottom plates 
and a Plexiglas tube with thickness 8 mm as the side- 
wall. The cell's inner diameter and height are both 18.5 
cm, so the aspect ratio of the cell is also unity. In the 
present work, 50-/im-diameter polyamid spheres (density 
1.03 gcm~ 3 ) were used as the seeding particles and a hor- 
izontal laser light-sheet with thickness ~ 2 mm was used 
to illuminate the particles in the mid-height plane of the 
cell. The measuring area is a square of 18.5 x 18.5 cm 2 , 
which covers completely the horizontal cross-section of 
the cell. The spatial resolution is 0.29 cm, corresponding 
to 63 x 63 velocity vectors in each 2D velocity map. De- 
note the laser-illuminated plane as the (x, ?/)-plane and 
the center of the mid-height horizontal cross-section of 
the cell as the origin O of the coordinates (see figure 
I13f a)). then two horizontal velocity components u(x,y) 
and v(x,y) are measured. In these PIV measurements, 
the cell was tilted by ~ 1° at position 1, and thus the 
orientation of the LSC was locked along the x-direction 
(see figure [lja)). For the PIV experiment, the measure- 
ments were made at Ra = 3.0 x 10 9 and 6.0 x 10 9 and 
at Pr = 4.3. A total of 20000 vector maps were acquired 
for each Ra at a sampling rate of ~ 2.2 Hz. As the 
two measurements give the same results, only results for 
Ra = 6.0 x 10 9 will be presented. 
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FIG. 11: Time-averaged veloc- 
ity vector map measured (a) 
near the top plate and (b) at the 
mid-height plane. For clarity, 
a coarse-grained vector map of 
size 32 x 32 is shown. The mag- 
nitude of the velocity \/w 2 + v 2 
is coded in both the color and 
the length of the arrows in units 
ofcm/s. Both time-averages are 
taken over a period of 151 min 
corresponding to 20000 velocity 
maps. 



To study the global motion of the central bulk fluid, we 
use the velocity vector V(t) — V x (t)x + V y (t)y to charac- 
terize the overall flow behavior at the mid-height plane, 
where V x (t) and V y (t) are the spatial averages of the local 
velocity components u(t) and v(t), respectively, within a 
circular area centered at the center of the measurement 
area. The diameter of the circular area is 5 cm and there 
are 221 vectors contained in this circular region for aver- 
aging. A circle with a diameter of 10 cm has also been 
used and the obtained results are essentially the same 
[3l| . The orientation 9 = arctan(V y /T4) and the magni- 
tude V = \V\ = \JV X + Vy of the averaged vector are 

then measures of the orientation and the strength of the 
motion of the central bulk fluid, respectively. 

Note that the horizontal velocity field investigated in 
Ref. [3l| was obtained near the top plate. Figure ITITa) 
shows the time-averaged velocity vector map measured 
at 2 cm from the top plate, which shows that the LSC 
is dominated mainly by the horizontal velocities. Thus, 
the spatial averaged velocity vector V(t) measured at the 
horizontal planes near the top plate can be used as a mea- 
sure of the orientation vector Vlsc of the LSC [3l| and 
its orientation 9 can be used to characterize the orienta- 
tion of the LSC at the corresponding height. However, 
the situation is different when the measurement is made 
near the cell's mid-height plane. Figure ITlT b) shows the 
time-averaged velocity vector map obtained at the mid- 
height horizontal plane of the cell. Because at the cell's 
mid-height plane the LSC is concentrated near the side- 
wall and dominated mainly by vertical velocities (see the 
circulation path of the LCS shown in figure QJb)), the ve- 
locity vector V(i) spatially averaged from the velocities 
in the central bulk region cannot be used to represent the 
orientation vector of the LSC and its orientation 9 is no 
longer a direct measure of the LSC's orientation. There- 
fore, 9 and </> m (the orientation of the LSC obtained from 
the azimuthal temperature profile) would exhibit differ- 
ent behaviors. As we shall see below, this is indeed the 
case. 

Figure [T2Ta) shows the PDF of the measured veloc- 
ity orientation 9. It is seen that the PDF of 9, which 



is significantly different from that of <fi m in figure [3tc), 
exhibits a bimodal distribution with two peaks located 
at the orientations (9 = ±0.25) that are perpendicular to 
the preferred orientation of the LSC (9 = 0). The proba- 
bilities of 9 — ±0.25 are nearly 3 times larger than those 
of 9 = and 9 = ±0.5, suggesting that the central bulk 
fluid is much more likely to move in the direction per- 
pendicular to the LSC's vertical circulation plane rather 
than in the LSC's preferred direction itself. To study the 
flow strength of the central bulk fluid in different orienta- 
tions, we calculate the conditional average (V\9) on the 
velocity orientation 9, as shown in figure [T2Tb). One sees 
that the averaged velocity magnitudes in the directions 
(9 = ±0.25) perpendicular to the LSC plane are much 
stronger than the magnitudes in all other directions, es- 
pecially than that in the preferred direction of the LSC. 
Figure [T2T c) shows the auto-correlation functions CV x 
(dashed line) and Cy (solid line). Both show oscillations, 
but the oscillation strength of V y (in the direction perpen- 
dicular to the LSC plane) is much stronger than that of 
V x (along the LSC's preferred direction) and V y oscillates 
more coherently than V x . In fact, a previous single-point 
velocity measurement has shown that at the cell center 
the strongest velocity oscillation is along the direction 
perpendicular to the LSC plane and that the strength of 
this oscillation decays away from the cell center towards 
the plates [l8l . [24| . We further note that the oscillation 
frequency of V y is the same as that of the off-center mo- 
tion of the LSC but 20% larger than that of V x . The 
existence of two different oscillation frequencies in tur- 
bulent convection with fixed control parameters (R a and 
Pr) have already been observed previously [H, 131(, but 
the reason(s) for the presence of two clocks in turbulent 
RB system remain unknown [33| . Here we offer a possible 
reason for the observed different oscillation frequencies 
of V x and V y . Even though the convection cell has been 
tilted by ~ 1° in the PIV measurement, the orientation 
of the LSC can still meander around the locked direction 
(the x-axis in the present case) over an azimuthal an- 
gular range of a few tens degrees [l|. The projection of 
the off-center oscillation onto the x-direction would thus 
generate the oscillation in V x , while the stochastic flue- 
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FIG. 12: (a) PDF of the orientation 9 of the spatially averaged 
vector V. (b) The conditional average (V\9) on the velocity 
orientation 9. (c) The auto-correlation functions of V x and 
V v . 



tuations of the projection angle would both reduce the 
coherence and change the frequency of the oscillation of 
V x . Note that the different oscillation frequencies for V x 
and V y have also been observed, but not explicitly recog- 
nized, by Refs. [18[ and [24| . Notwithstanding the above, 
the overall picture emerging from the PIV measurements 
is that the translational motion of the central bulk fluid 
in the direction perpendicular to the LSC's vertical cir- 
culation plane is more probable, much stronger and more 
coherent than the motions in all other directions, which 
provides a direct evidence for the off-center motion of the 
LSC at the mid-height plane. 

In a previous work, Sun et al. [24] measured the 2D 
velocity field in the (y, z)-plane, which is perpendicular 
to the LSC's vertical circulation plane and studied the 
phase relationships among the horizontal velocity com- 
ponents v(y) along the y-axis and among v(z) along the 
z-axis. For v(y) along the y-axis, their results showed 
that the velocity component v(y) at different values of 
y are highly correlated and have a common phase across 
the cell's entire diameter. For v(z) along the z-axis, their 
results revealed that the horizontal velocity component 
v(z) at different height z remain in phase along the z- 
axis mainly in the middle one-half of the cell, while the 
horizontal velocity components v(z) near the upper and 
lower conducting plates gradually lag behind those in the 
central region of the cell. Both of these results imply that 
the motion of the central bulk fluid is spatially coherent 
and the off-center oscillation is the dominant mode of the 
central bulk fluid motion. 

The bulk off-center mode of the flow field can also 
be revealed by investigating the variation along the x- 
axis of the y velocity component v(x, 0). Here we study 
the relationship between v(x, 0) and v(— x, 0) for vari- 
ous values of x(> 0). If the azimuthal rotation is the 
dominant motion of the LSC, v(x 1 0) should anticorre- 
late with v{— x,0), i.e., v(x,0) is along the y-direction 
when v(—x,0) is along the — y-direction and vice versa. 
While v(x, 0) correlates strongly with v(— x, 0) for the sit- 
uation that the off-center motion is the LSC's dominant 
mode, i.e., v(x, 0) and v(— x, 0) are along the same direc- 
tion. Figure [T37b) shows the typical time series of v(x, 0) 
and v(—x,Q) for x/D = 0.11. It is seen that v(x,0) and 
v(— x, 0) are very similar to each other and along the same 
direction for most of the time. In addition, both v(x, 0) 
and v(—x, 0) are found to exhibit a well-defined periodic 
oscillation, corresponding to the off-center oscillation of 
the LSC. To characterize quantitatively these features, 
figurefTlTa) shows the cross-correlation functions C(r, x) 
between v(x,0) and v(— x,0) for x/D = 0.11 (solid line) 
and 0.33 (dashed line). It shows that both functions os- 
cillate and have a large positive peak located at t = 0, 
indicating that v(x,0) correlates strongly with v(—x,0). 
It is further found that this peak exists for all x and 
all horizontal velocity oscillations of v(x,Q) have a com- 
mon phase across the cell's entire diameter. Figure [14T b) 
shows the cross-correlation coefficient C(r = 0, x) as a 
function of x. One sees that C(r — Q,x) > 0.2 for all x 
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FIG. 13: (a) Coordinate system and 
schematic diagram of the mid-height hor- 
izontal cross-section of the cell. The ori- 
gin O of the coordinates is chosen as the 
center of the cross-section and thus the z- 
axis (not shown in the figure) is along the 
cell's central vertical axis. The two verti- 
cal arrows show the velocity components 
v located at positions (a;,0) and (— x,0) 
and the dashed line shows the ir-axis. (b) 
Time traces of v(x, 0) and v(— x, 0) for 
x/D = 0.11. 
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0, x) as a function of x. 



and it has a peak near the sidewall (around x/D ~ 0.42). 
The peak of C(r = 0, x) near the sidewall represents the 
anticorrelation between w/ ljm and w Cjm , i.e. the off-center 
motion of the LSC, as illustrated in figures [TWb) and (c). 
Both the strong positive peak located at r = and the 
oscillation of the cross-correlation functions C(r, x) be- 
tween v(x, 0) and v(x, 0) for all x reveal the fact that the 
dominant motion of the LSC at the mid-height plane is 



indeed the off-center oscillation, which provides another 
direct evidence for the off-center mode of the LSC at the 
mid-height plane. 



IV. CONCLUSION 

To conclude, we have made a detailed experimental 
study of the oscillations of the large-scale circulation 
(LSC) in a cylindrical turbulent Rayleigh-Benard convec- 
tion cell with aspect ratio unity using water as working 
fluid. Direct spatial measurements of both the temper- 
ature and velocity fields were carried out to study the 
motion of the LSC. 

Direct measurements of the horizontal velocity field at 
the cell's mid-height plane using the PIV technique were 
made at Ra = 3.0 x 10 9 and 6.0 x 10 9 and Pr = 4.3. It 
is found that the horizontal translational motion of the 
central bulk fluid in the direction perpendicular to the 
vertical circulation plane of the LSC is more probable, 
much stronger and more coherent than the motions in 
all other directions. For the velocity component v(x, 0) 
perpendicular to the LSC plane, v(x, 0) for all values of x 
are found to correlate strongly with each other, especially 
for the velocities near the sidewall, and have a common 
phase across the cell's entire diameter. Both of these 
results provide direct evidences for the off-center mode 
of the LSC at the mid-height plane. 

The temperature measurement was performed over the 
Rayleigh-number range 9 x 10 8 < Ra < 6 x 10 9 and at 
fixed Prandtl number Pr — 5.3. At each Ra the instan- 
taneous azimuthal temperature profile along the sidewall 
at three different heights, i.e. H/A, H/2, and 3H/4 from 
the bottom plate, were measured. To study the motion of 
the LSC, we developed a new method, the temperature- 
extremum-extraction (TEE) method. Using this method, 
the azimuthal angular positions of the hot ascending and 
cold descending flows at each height were obtained by 
making a quadratic fit around the highest and lowest 
temperature readings, respectively, and the line connect- 
ing these two positions is the central line of the LSC 
band. The orientation of this line is thus the orientation 
of the LSC and the distance between this line and the 
cell's central vertical axis is defined as the off-center dis- 
tance. The motion of the LSC is therefore decomposed 
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into two different modes: the azimuthal mode and the 
translational or off-center mode. When compared to the 
sinusoidal-fitting (SF) method, it is found that, as far as 
the azimuthal motion of the LSC is concerned, the two 
methods give the same information in terms of both the 
orientation and the flow strength of the LSC, while the 
TEE method can in addition determine the LSC's off- 
center motion that is missed by the SF method. With 
the TEE method, the LSC's central line is found to os- 
cillate time-periodically around the cell's central vertical 
axis with an amplitude being nearly independent of the 
turbulent intensity, leading to the off-center oscillation of 
the LSC. 

Results obtained using the TEE method further reveal 
that both the torsional and off-center modes of the LSC 



are manifestations of the out-of-phase azimuthal posi- 
tional oscillations of the hot ascending and cold descend- 
ing flows. When these two flows are near the top and 
bottom plates, the resulting motion is torsional. When 
they are near the mid-height plane, the resulting motion 
is the off-center oscillation. 
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